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Abstract 

Thomson, D. J. & Devenish, B. J. [J. Fluid Mech. 526, 277 (2005)] and others have 
suggested that sweeping effects make Lagrangian properties in Kinematic Simulations 
(KS), Fung et al [Fung J. C. H., Hunt J. C. R., Malik N. A. & Perkins R. J. J. Fluid 
Mech. 236, 281 (1992)], unreliable. Here it is shown through a novel analysis based 
upon analysing pairs of particle trajectories in a frame of reference moving with the 
large energy containing scales of motion that the normalized integrated error e T K in the 
turbulent pair diffusivity (AT) due to the sweeping effect decreases with increasing pair 
separation (cq), such that e —> 0 as cq/q —> oo; and e J K —> oo as cq/q —> 0. q is the 
Kolmogorov turbulence microscale. There is an intermediate range of separations 
1 < cq/q < oo in which the error e J K remains negligible. Simulations using KS shows 
that in the swept frame of reference, this intermediate range is large covering almost the 
entire inertial subrange simulated, 1 < eq/q < 10 5 , implying that the deviation from 
locality observed in KS therefore cannot be atributed to sweeping errors and could be 
real. This is important for pair diffusion theory and modeling. 


PACS numbers: 47.27.E?, 47.27.Gs, 47.27,jv, 47.27.Ak, 47.27.tb, 47.27.eb, 47.11.-j 

Keywords: Turbulence, diffusion, particle pair, pair diffusivity, Kinematics Simulations, 
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Introduction 

Turbulent particle pair diffusion has attained somewhat of an iconic status in the 
turbulence community, many researchers having addressed this topic over the decades. 
Nevertheless, most if not all theories of turbulent particle pair diffusion in homogeneous 
turbulence with extended inertial ranges have been based upon the hypothesis of 
locality since Richardson in 1926 l], and Obukov in 1941 [2], 

Richardson pioneered this field and introduced the idea of a scale dependent pair 
diffusivity as the fundamental quantity of interest in turbulent pair diffusion studies. 
The turbulent pair diffusivity is defined as, 
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where 1 (f) is the pair displacement vector at time t, l = | 1 |, v(Z) is the pair relative 
velocity, and (■) is the ensemble average over all particle pairs. 

The locality hypothesis can readily be applied to generalized power law spectra of the 

type, E(k ) ~ k p , for 1 < p < 3; the pair diffusivity then scales like K(l,p) ~ a l p with 
7 p = (1 +p)/ 2, 31, where a 2 = (l 2 ). For Kolmogorov turbulence p = 5/3, this gives the 

well known Richardson scaling K ~ er 4 ^ 3 , which is equivalent to ( l 2 ) ~ f 3 j2j. 

Kinematic Simulations (KS) [4} 5] has often been used to investigate turbulent pair 
diffusion, even though it does not yield the assumed locality scaling; for p = 5/3, KS 
gives 7 Koi ~ 1-53 > 4/3. For this reason, it has been widely assumed that KS must be in 
error. Thomson & Devenish |6 argue that the turbulent pair diffusivity must scale like, 

Km) ~ s(i)T„m), ( 2 ) 

where S(l) is the structure function of the turbulence velocity field and t s (1) is an 
effective time scale of velocity increments. In real turbulence, assuming locality scaling 
for S(l) ~ Z 2 / 3 and for t s (1) ~ l 2 / 3 leads to Richardson’s classical scaling K(l) ~ ct ; 4 ^ 3 , 
where we evaluate K at typical values of /, namely cq = (l 2 ) 1 / 2 which is commonly 
assumed in these studies. 

Thomson & Devenish argue that in KS because of the lack of true dynamical sweeping, 
the time scale must be a sweeping time scale, namely r s ~ l/U s which is a time scale for 
the large scales flow to cut through smaller local eddies, U s being the sweeping velocity 
scale. This leads to, K ~ (t / 5/3 . Even when the rms turbulence velocity u' is taken 
instead of U s , they obtained K ~ cr 44//9 . 

They concluded that whereas locality is true in real turbulence, it is not true in KS. In 
turbulence the large energy containing eddies carry the smaller eddies, but in KS as 
there is an absence of true dynamics the large scales force the fluid particles to cut 
through the smaller eddies in an unphysical manner, a view supported by Nicolleau & 
Nowakowski 17 , and Eyink & Benveniste [8]. 

However, Thomson & Devenish’s scaling argument leading to equation (2) adresses only 
the scaling laws for K, but does not quantified the actual errors in the diffusivity K in 
KS - is it large or small ? It is prudent, therefore, to seek an alternative, a more 
analytic, approach to address this question, which is the main concern of this work. 

Here we re-examine the sweeping effect in KS with a view of quantitfying the error in 
the KS pair diffusivity K s compared to the physical pair diffusivity K. For this purpose, 
we focus upon the differences in the relative velocities along pairs of particle paths in the 
sweeping frame of reference. This frame of reference accounts for the physical sweeping 
effect of the largest energy containing scales; but a residual sweeping effect still remains 
due to the largest inertial range eddies sweeping the smaller inertial range eddies. 

Consider Fig. 1 which shows a particle pair with separation 1 in the inertial subrange 
being swept by a large scale flow. A real fluid particle pair will be swept by the physical 
velocity field u and will follow certain particle paths; but a KS flow will transport the 
pair along neighbouring particle paths due to an additional KS sweeping motion, u s , 
and thereby force the particles to cut through local flow structures. 

The large scale physical sweeping velocity are assumed not to affect the relative motion 
of particles in a pair in the inertial subrange. The critical question is, are the deviations 
from the physical trajectories induced by KS in the pair diffusion process large or small? 

In the ensuing analysis, it is the error between the KS and the physical pair diffusivities, 
| K s — K |, that will be calculated. We will consider generalised power law energy 
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spectra, E(k) ~ k~ p , because the analysis must be valid for all such power spectra and 
this will add weight to the results and conclusions that can be drawn from this work if 
validated over the wholw range of p considered. 

The main questions of interest are, is the KS sweeping error large or small, and in what 
range of separations? These questions are addressed first through a novel mathematical 
analysis focussing upon pairs of neighbouring particle trajectories. This is then verified 
against simulations using KS with very large inertial subranges. 

In Section 1, we derive an expression for the error in the pair diffusivity in KS flows by 
analysing neighbouring trajectories in the swept frame of reference. In Section 2, the KS 
method is discussed and simulation results presented. In the final Section 3, we discuss 
the results and its implications for theory and modeling. 


1 The normalized error in the pair diffusivity 

1.1 The numerical timestep error 

In the swept frame of reference, the relative motion of a particle pair in the inertial 
subrange is unaffected by the sweeping action itself. This can be mimicked in KS by 
setting E{k) =0 for k <k\. However, there still remains a residual sweeping caused by 
the largest of the inertial scales sweeping the scales local to the pair separation. 

Consider an ensemble of particle pairs released in a field of homogeneous turbulence at 
time t = 0 with some small initial separation 1$. At some time t later, the ensemble 
average of the separation is assumed to be well inside the inertial subrange and the 
relative motions are independent of 1 $ j9]. 

Consider the particles in one of these pairs, labeled 1 and 2, as shown in Fig. 1. The 
particle locations are xi(t) and X 2 (t) respectively at time f; and the pair displacement is 
1(f) = X 2 — xi, and l{t) = |x 2 — xi|. The physical flow is u(x,t). All quantities are 
assumed at time t unless otherwise stated. 

At time t the additional (or residual) KS sweeping flow u s (x,t) is ’swiched on’ - this is 
not to be confused with the total KS velocity field which is (u + u s )(x,t), see Fig. 1. 

The physical flow u(x, t) transports the particles to xi(f*) and X 2 (t*) respectively at 
the next time step t* = t + dt: while the KS flow (u + u s )(x, t) transports the particles 
to x*(t*) and x^t*) respectively. Note that I s = Xg — xf, and I s = |x| — x®|. 

The superscript * will refer to quantities at time t*, e.g. I* = l(t + dt). The superscript 
s will refer to quantities related to the KS residual sweeping, e.g. I s (t*) = I s (t + dt). 

The following quantities are defined: 

u = u(x, t) is the physical fluid velocity field 

u s = u s (x,f) is the additional (residual) sweeping velocity field 

v(l) = u(x 2 ) — u(xi) is the physical relative velocity 

v s (1) = u s (x 2 ) — u s (x!) is the additional (residual) relative velocity 

u = (u + u s )(x,t) is the total KS velocity 

v(l s ) = v(l s ) + v s (l s ) is the total KS relative velocity 

Simplifying the notation as much as possible, e.g. u 2 = u(x 2 ,t), and Uj = u(x 2 ,t + dt), 
yields 
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Figure 1 . Schematic diagram illustrating the system discussed in the text. The 
locations of two nearby particles, labelled 1 and 2, are located at Xi(f) and x 2 (t) 
respectively, with turbulence velocities Ui(f) and u 2 (f), at time f; their separation is 
l(t) = |x 2 — xi|. They are transported with velocities (ui + u|)(f) and (u 2 + u|)(f) 
respectively to the new locations x s (f) and x|(t) at the next time step t + dt, as shown. 


* 5(0 

= *i{t*) + ul(t)dt 

(3) 

* 5(0 

= x 2 (t*) + U-2(t)dt 

(4) 

JS* 

= x£(t*)-xf(f) 



= 1 (f) + (u| — u {)dt 



= l*+v s (l )dt 

(5) 


v(Z s *) is calculated at the new KS swept particle locations. Using Taylor expansions 
wherever necessary, assuming that the velocity fields are at least twice differentiable in 
space and at least once in time, 


v(n = (u + u s )(x^*))-(u + u s )(x s (;£*)) 

= u(x*(f*)) - u(x s (U)) + u s (x s (f*)) - u s (x s (U)) 

= v(T) + v s (1*) + 

(u^ • Vu 2 (f*) — uj • Vui(t*)) dt + 

(uj • Vu|(t*) — uj • Vu {(t*)) dt 

+0(dt 2 ) ( 6 ) 

The pair diffusivity at time t* is, K * = (1* ■ v(Z*)} - we ignore constants of 
proportionality, like 2, because we are interested only in the power scalings in this work. 
The KS equivalent is K s * = (I s * • v(Z s *)). Using equations (5) and (6) and ignoring 
terms of order dt 2 and higher, 
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K s * » (1* • v(l*)) + (1* • v s (l*)) + 

(1* • (u s 2 (t) • Vu 2 (f) - u f(t) • V Ul (t*)))dt + 
(1* • (u5C t) ■ Vu5(f) - uf(t) • Vu;(t*)))dt + 
(v a (i)-v(r))dt + (v»(i)-v 8 (r))dt 


The timestep error between the KS and physical diffusivities for a given timestep dt is, 
E k = | K s * — K*\. Using the expansion u s (x 2 (f)) « uf + 1 ■ Vuj in equation (7), yields 


e k « (r • v s (i*)) + 

(r • ( u ; • v)v(i*))dt + 

(T • (uj ■ V)v 8 (l*))dt + 

(1* • (1- VK- Vu 2 (t*))dt + 

(1* • (1-VK-Vu 8 (t*))di + 

(v 8 (i)-v(r))dt + 

(v 8 (l) ■ v 8 (l*))dt (8) 

The time scale of the sweeping T s is much larger than the local time scale of the pair 
separation. Hence, in the last four terms 1(f) is replaced by l(f*) without affecting their 
magnitudes or scalings (the associated errors are ~ 0(dt 2 ) which is neglected). 

All the terms in equation (8) are now evaluated at the same time f*, so without loss of 
generality t* is replaced by t and the superscritp ’*’ is dropped. The subscript V is also 
dropped because of homogeneity. Equation (8) now simplifies to, 


E k « <l-v s ) + 

(1- (u s ■ V)v)df + 

(1 • (u s • V)v s )dt + 

(1 • (1 • V)u s • Vu 2 )dt + 

(1 • (1 • V)u s • Vu 2 )dt + 

(v s • v)df + 

((v s ) 2 )dt (9) 

It is reasonable to assume that the (residual) sweeping flow field u s , which is caused by 
the largest of the inertial range eddies, is close to uniform across small distances, and 
therefore the relative velocities across local scales l that it induces is small. However, 
gradients of the relative velocity v(l) itself can be large. The magnitude u s = |u s | is 
assumed large compared to v(l) = |v(Z)|, and also as compared to v s (l) = |v s (Z)|. u a 
scales differently to v(l). 

v(l) also scales differently to v s (l), the former being governed by inertial range 
turbulence scaling, and the latter by differences in the residual sweeping velocity across 
a small distance l. This can be seen clearly in the limit of uniform (parallel) sweeping 
flow, where the v(Z) is unaffected, but v s (f) = 0 and all the terms on the right hand 
side in equation (9) are zero except for the second term. This indicates that the second 
term in equation (9) makes the dominant contribution to the error. 
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Consider generalized energy spectra of the form E(k) = e 2 ^ 3 L 5 ^ 3 ~ p k~ p , for ki < k < k v 
and for 1 < p < 3, and with k v /ki 1. In the swept frame of reference, E(k) = 0 for 
k < k\. The rate of energy dissipation is e ~ U 3 /L , where U is the velocity scale in the 
energy containing scales. 

The previous discussion implies that |Vv s | -C and therefore, 

v a (l ) « (£) U. (10) 

The energy in turbulent inertial scales local to l is, v 2 {l) ~ E(l/l)/l, and therefore, 

p —1 

v(l) ~ (l) 2 U (11) 

p —3 

and, |Vv(0| ~ (jr'j j (12) 

It is usual to assume the scaling l ~ cr; as previously mentioned. Then, the second term 
in equation (9) is given by, 


p— 1 

E 2 = (1 • (u s ■ V)v)dt ~ 2 Uu s (l)dt , (13) 

All the other terms in equation (9), labeled respectively E\, E 3 ,..., E^, scale 
proportional to u s or are much smaller, and this leads to the following estimates, 


Ei 

E 2 

e 3 

e 2 

Ei 

E 2 

E5 

E 2 

Eg 

E 2 

E 1 

E 2 



(14) 


All of these ratios are small for <JilL < 1 and 1 < p < 3. This is true even in the 
expression for E- t /E 2 because the factor L/Udt is subdued by the very small factor in 
the brackets. It is reasonable to conclude that the 2nd term in equation (9) is dominant 
and therefore Ek « E 2 . 

To estimate E 2 itself, an estimate for u s (l ) is needed. The inertial subrange contains 
only about 1% of the the total energy in the turbulence across the entire wavenumber 
range, 0 < k < oo, such as in a von Karman spectrum E v k , i.e. E^ s ss E v k/ 100. This 
means that the root mean square turbulent velocity in the inertial range is 
approximately, Uks ~ U/ 10, where E v k ~ U 2 . 
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x 

Figure 2. The function f{x) = y/l — x p ~ 1 x ( - p ~ 1 ^ 2 , 0 < x < 1, for selected powers 
1 < p < 3. 

The wavenumbers that contribute to the sweeping of particle pairs at separation cq are 
in the range k\ < k < ki, where ki ~ l/cq. In the KS sweeping frame of reference, the 
larger inertial scales are the sweeping scales, and the energy in these scales is 
approximately, 


(u s ) 2 ~ [ 1 E{k)dk = [ 1 a k e 2/3 L 5 / 3 - p k- p dk (15) 

Jki Jki 

and using e ~ ( Uk s ) 3 /L ~ (U/10) 3 /L, 


Using this in equation (), the residual error between the physical and the KS pair 
diffusivities in the inertial subrange of pair separations in the swept frame of reference, 
per unit timestep (retaining Ek to represent this quantity), is 

B K W~E, = (17) 

Cfc is the constant of proportionality, which can depend upon p. 

For p = 5/3, the residual error per unit timestep is, 


this gives 


U 

10 


1 - 


(?) 


P-1 


(16) 
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(18) 


E k ( 5/3) 



1/3 U 2 

10 ' 


As p —► 1, E# —> Eif(l) ss CkU 2 /L ss Constant for ai/L < 1, but is nearly zero close 
to <Ji/L = 1. In this limit, p —> 1, the pair diffusion is strongly local and is not affected 
by long range sweeping. 

For p = 3, the residual error per unit timestep is negligibly small for cp/E < 1, 


Ek{ 3) « ^«^ ( 5 /3) - (19) 

In this limit, nearly all the energy is contained in the largest scales and inertial 
subrange scaling is no longer applicable. 


1.2 The integrated error 

Equations (17) provides a way of estimating an upper bound for the integrated residual 
error. Fig. 2 shows the log-log plots of the factor f(x) = \/\ — x p ~ 1 x^ p ~ 1 ^ 2 for selected 
powers in the range 0 < p < 3. The range over which \/l — x p_1 ss 0 is very short and 
close to x = 1; but over the rest of the range \/l — x p -1 « 1. Hence, to a good 
approximation, f(x) < x (p _ 1 V 2 f or a u x < 1, and using this in equation (17) with 
x = ai/L yields, 


Ek(p) < ■ (20) 

The integrated residual error, E^-, over a period of time is, 



Assuming the pair separation scaling a 2 ~ t Xp , where t is the time and for some Xp > 0j 
yields, 


E t k < C k 


UL 

TiT 



( 22 ) 


If the pair diffusivity scales like, I\ ~ a' / fp , for some > 0 then \ p — 2/(2 — y p ) is an 
exact relation. 

The most important quantity is the normalised integrated residual error with respect to 
the pair diffusivity, e£-. Using the above expression for \ pi and replacing the scalig with 
L by scaling with ip and absorping all ensuing constants in to a new constant Ak, leads 
to 
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c / -Ejl< 

e K- K < 


10[ - 


27p - 3 


(23) 


For strict locality scaling = (1 +p)/ 2, and this becomes 


-K 


< 


Ah, 


10 


7 P -1 ‘ 


(24) 


Since y p — 1 >0, e J K decreases with increasing pair separation for all p > 1. For 
p = 5/3, we have 7 P = 4/3 and we obtain, 





(25) 


Thus, as (Ji/rj —> oo then e ! K decreases; and as cq/r? —> 0 then increases. 

Even for non-local scaling, assuming that 7 P does not deviate too far from the local 
scaling, the above order of magnitude for e ! K is still approximately true. In fact, in KS 
we know, Section 2.3, that 7 P is slightly greater than locality so the errors will be 
slightly smaller than in eqution (25). 

The magnitude of these errors will also depend upon Ak\ but Ak cannot be determined 
from theoretical considerations alone, although it is reasonable to assume that Ak is 
small enough for the e ! K to be negligible towards larger separations as ai/rj —>• oo - this 
assumption will be justified by numerical simulations. Even so this does not guarantee 
that the error will remain small at smaller separations, since e^- —> oo as cq /77 —> 0. 

However, between these two asymptotic limits there must exist an intermediate range of 
separations, between 1 < ai/p < oo, where the errors remain negligible. The crucial 
question is, how wide is this range of scales? 

If this intermediate range is short then the KS errors will be significant at almost all 
separations. However, if it is long then the KS errors will be negligble inside this 
intermediate range where true inertial subrange scaling can be expected to be 
manifested. 

To determine the extent of this intermediate range of scales, if it exists at all, 
simulations with KS must be performed with very large inertial subranges. 


2 Simulations and Results 

The normalised integrated error, e T K , is scale dependent and reduces with increasing 
separation. The KS diffusivity is given by, K s ss K( 1 + e^) —> K as crj/ry —>■ oo. It is 
expected that if there is an appreciable intermediate range where the errors are 
negligible, then the power scaling in K s must be constant and asymptotic to the 
limiting case where ai/p -4 oo. The extent of this intermediate range is determined by 
the range over which the power scaling in K s is constant. 
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Furthermore, significant levels of the sweeping error means that the fluid particles cut 
through KS eddies, and therefore must be accompied by high levels of noise - the larger 
the relative sweeping error the larger the noise level. 

Thus, where the errors are negligible it is expected that the correct power law scaling, 
K s k, K, will be observed in that part of the of the inertial subrange. 

On the other hand, where the errors are significant it is expected that K s will deviate 
from the true power law scaling for K and also be accompanied with significant 
statistical noise due to fluid particles being swept through local eddies in that part of 
the inertial subrange. Even in this case, however, it is expected that the errors and the 
associated noise diminish as the pair separation increases. 

2.1 Frames of reference 

Comparison will be made between two cases: first, where K s is obtained from KS in the 
physically correct sweeping frame of reference; and second, the case where large scale 
random sweeping velocities are explicitly added to the flow. 

The analysis for the latter case is similar to that which leads to equation (24), except 
that the factor of 10 in the denominator is now 1. Both of these cases can therefore be 
written collectively as, 


e 


i 

K 



(26) 


where C = 10 in the swept frame of reference, and C = 1 when large random scales are 
included in the simulations - the residual errors are an order of magnitude smaller in 
the swept frame of reference. 

Case 1: Swept frame. Set the spectrum to be E{k ) ~ k~ p in the inertial subrange, 
and set E(k ) = 0 for k < k\. 

Case 2: Non-swept frame. Set the spectrum to be E{k) ~ k~ p in the inertial 
subrange as in Case 1, and add E(k) = Eo5(k — ko ) at some low wavenumber ko < k\, 
and such that Eq is the energy in the von Karman spectrum in the range 0 < k < k\. 

A very small fixed timestep, smaller than any timescale in the system dt -C t v , is used 
in all the simulations reported here. t v is the Kolmogorov time scale. 


2.2 Kinematic Simulations 


Kinematic Simulation [4,5] is a Lagrangian method for particle diffusion in which the 
velocity fileld is prescribed as a sum of energy-weighted Fourier modes. It is akin to the 
widely used random flight type of statistical models in which the dynamical interactions 
between turbulent length scales is not explcitly modeled, rather the overall effect on the 
statistical moments of particle diffusion is mimicked. In KS this is accomplished by 
specifying the energy spectrum E{k). KS continues to be used in turbulent diffusion 
studies for both passive and inertial particle motion, including cases with generalized 
power-law energy spectra of the form E(k) 

Malik & Vassilicos 


10 , Turfus 


11 


Fung & Vassilicos 12 


k p for p > 1, Maxey 
[13]. Meneguz & Reeks 14T~carried out a 


DNS of inertial particle motion, and compared it to results from KS which they found 
to agree well with the DNS. 
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<V n 

Figure 3. The turbulent diffusivity as log (K/(r)V v )) against log (cr;/r?) obtained from 
KS. From top to bottom, p = 1.01, 5/3, 3. Case 1 (red lines), Case 2 (green lines). 
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KS generates turbulent-like non-Markovian particle trajectories by releasing particles in 
flow fields that are incompressible by construction and which satisfy Eulerian statistics 
up to second order. A turbulent flow field realization is produced as a truncated Fourier 
series, 


N 

W (x, t) = [(A n x k n ) cos (k n • X + co n t) + (B n x k n ) sin (k n ■ x + uj n t) 

n—1 


(27) 


where N is a suitable number of representative wavemodes, typically hundreds for very 
long spectral ranges, k v /ki 1 . k„ is a random unit vector (k„ = k n k n and 
k n = |k„|). The coefficients A„ and B„ are chosen such that their orientations are 
randomly distributed and uncorrelated with any other Fourier coefficient or 
wavenumber, and their amplitudes are determined by (A 2 ) = (B 2 ) oc k n E{k n ), where 
E(k), k\ < k < k v , is the turbulent energy specturm. The angled brackets (•) denotes 
the ensemble average over many flow fields. This construction ensures incompressibility 
in each flow realization, V • u = 0. The flow field ensemble generated in this manner is 
statistically homogeneous, isotropic, and stationary. 

An important feature of KS is that unlike some other Lagrangian methods, by 
generating entire kinematic flow fields in which particles are tracked it does not suffer 
from the crossing-trajectories error which is caused when two fluid particles occupy the 
same location at the same time in violation of incompressibility; but because KS flow 
fields are incompressible by construction this error is completely eliminated. 

The energy spectrum E(k) can be chosen freely within a finite range of scales. In 
turbulent particle pair studies the interest is in Kolmogorov-like power law spectra, 


E(k) = C E s 2 / 3 L 5 ^- p k- p , h<k< k v {= 2ir/ri), 1 < p < 3 (28) 

C E is a constant. The largest represented scale of turbulence is 27 r/fci, and the smallest 
is the Kolmogorov scale r) = 2ir/k v . The constant is normalized such that the total 
energy contained in the range k\ < k < k^ is 3(u , ) 2 /2, where u' is the rms turbulent 
velocity fluctuations in each direction. e{p) is determined by integrating the spectrum, 
jkr, = 3(u') 2 /2. (p = 1 is a singular limit which is not consider here.) 

v v = ( er z ) 1 / 3 is the velocity micrcoscale, and t v = e -1 / 3 ?^ 2 / 3 is the Kolmogorov time 
micrcoscale. 

The frequencies are chosen according to usual practice to be proportional to the 
eddy-turnover frequencies, i.e. uj n = A \/k^E(k n ). The choice of A is somewhat 
arbitrary, but provided A < 1 it does not affect the diffusion scaling itself - even frozen 
field with A = 0 yields the same scaling [15]. A = 0.5 is a common practice in KS which 
is also chosen here. 

The distrbution of the wavemodes is geometric, k n = k\r n ~ 1 , with r = (fc^/fci) 1 ^^ -1 ^. 
The grid size in wavemode-space of the n th wavemode is Sk n = k n (y/r — l/y^). 

A particle trajectory is obtained by integrating the Lagrangian velocity Wz,(t), 

(}-v 

_= Wi (i) = W(x,i). (29) 

Pairs of trajectories are harvested over a large ensemble of flow realizations and pair 
statistics are then obtained from it for analysis. 
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The turbulent difusivity itself can be computed in two ways. Directly from the formula 
K(l) ~ (1 ■ v(Z)), i.e. the ensemble average of the scalar producted of v and 1 . But it has 
been found that using the equivalent formula, K(l) ~ d(l 2 )/dt , i.e. the derivative of the 
(l 2 ), converges faster statistically needing a much smaller ensemble of trajectories, 
although the two methods give identical results for large enesmbles of particle 
trajectories. The latter method has been adopted here. 

Lagrangian statistics are the physically meaningful output from KS. It is not correct to 
compare the kinematically generated flow fields directly with DNS flow fields. As such, 
KS is like Lagrangian methods such as Random Walk models where an individual 
particle trajectory has no physical meaning, but the ensemble average over many such 
random trajectories produces physically meaningful Lagrangian statistics. 


2.3 Results 

KS simulations were performed with L = 1, k± = 1/L = 1, and k v = 10 6 , Ce = 1.5 
(Kolmogorov constant) and u' = 1 . There were 200 wavemodes per realization. 

In Case 1 (swept frame of reference), E(k) = 0 for k < 1. 

In Case 2, large scale random sweeping were added at the low wavenumber kg = 1/10, 
with E(k) = E 0 5(k — k 0 )- The energy in these sweeping scales, E 0 , was equal to the 
energy contained in the von Karman turbulence spectrum for k < 1. ko corresponds 
approximately to the location of the peak in the von Karman spectrum. 

In both cases, three different power spectra were considered with, respectively, 
p = 1.01,5/3, and 3. With 8 pairs released in 5000 flow realizations, the Lagrangian 
statistics were obtained from 40,000 particle pair trajectories. 

Fig. 3 shows log-log plots of the pair diffusivity K/(rjv v ) against eq/77, for Case 1 (red 
lines) and Case 2 (green lines). The energies in the two cases are different, so Case 2 
plots have been shifted vertically in order to compare the two cases directly. This does 
not affect the scalings (the slopes) which is the main interest here. Hence the ordinate is 
shown without scale. 

For p = 1.01, the two cases align with a constant power-law scaling, 71.01 ~ 1.07, over 
most of the inertial subrange of scales and there is very little statistical noise, indicating 
that e I K (ai) <C 1 at all separations in this part of the inertial subrange in both cases. In 
this limit, locality is very strong, and the relative motion is unaffected by the long range 
sweeping. The obtained slope is indeed very close to the exact locality scaling of 1.005. 

For p = 5/3 (Kolmogorov turbulence), in Case 1 (red) a clear power-law scaling is 
observed, 7 Koi ~ 1-53 > 4/3, and very little statistical noise in the range 
1 < ai/r] < 10 5 , indicating that e I K (ai) -C 1 in this range of scales. Case 2 (green) 
deviates increasingly from Case 1 at small inertial separations where it is also 
accompanied with increasing levels of noise. Nevertheless, the agreement between the 
two cases for cq/77 > 10 2 is good. 

For p = 3, the two cases overlap with a power-law scaling, 73 = 2, with almost no 
statistical noise. In this limit nearly all the energy is in the large scales and inertial 
range scaling is no longer applicable; rather uniform strained motion with the 
characteristic slope of 2 is obtained. 
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3 Discussion and Conclusions 


All the results in Fig. 3 are consistent with the numerical analysis and the theoretical 
predictions in section 2. The constant power law scaling over most of the inertial 
subrange of separations, 1 < <ri/r] < 10 5 , and the very low level of statistical noise in the 
swept frame of reference are especially important. (The departure from this for ai/p < 1 
observed in Fig. 3 is outside of the inertial subrange.) 

The KS sweeping errrors in this frame of reference is thererfore negligible for most 
practical purposes. It is possible that KS could produce negligible sweeping errors in 
even bigger intermediate ranges than reported here, but the current simulations are the 
maximum size of inertial subrange possible, k n /ki = 10 6 , with double-precision 
accuracy. 


It is remarkable that even when large scale sweeping is included, the KS sweeping errors 
remain small in the range cri/rj > 10 2 . 


It is also noted that some Direct Numerical Simulation (DNS) show pair diffusion which 
appear to display locality scaling, see 16 for example. However, the maximum inertial 
range obtained in DNS to date is around k v /k\ ~ 10 2 , which is much shorter than 
required to test pair diffusion scaling reliably - that requires k v /ki > 10 4 . Thus the 
current KS results cannot be compared directly with DNS at the present time. However, 
for low Reynolds KS has been validated against DNS for turbulent pair diffusion by 
Malik & Vassilicos 13 ; here not only did the pair diffusion from KS closely match the 
DNS results with the same energy specturm, but the fourth order statistic, the kurtosis 
in the pair separation, also matched remarkably well. 


The main contribution of this work is that it has been demonstrated that in a reference 
frame moving with the large energy containing scales the sweeping errors in the 
turbulent pair diffusion process in KS is negligible in the inertial subrange where, 

1 < ai/rj < 10 5 . This is significant not only because it amends the previous theory 
of [6]j8] , but it is important for turbulent diffusion modeling and applications in 
general, and especially for pair diffusion studies. 

If the sweeping errors are negligible, why is locality scaling not oberved for p = 5/3 
where KS yields 7 Koi ~ 1-53 > 4/3? There are two possible explanations. First, other 
errors could be present in KS, not due to the sweeping, but to as yet unknown sources; 
but no one has proposed any such source of error in, and so this remains very 
speculative. 


The second possibility is that the locality hypothesis itself may be error. This goes 
against the widely accepted theory of Richardson [1], and no one has proposed 
alternative theory. On the other hand, it should be noted that the locality scaling for 
the pair diffusion has never been confirmed unequivocally as noted by Salazar & 
Collins 


17 so there is room for new thinking in this field. 


These are currently the subjects of active research by the author. What can be said for 
certain at the present time is that the departure from locality scaling observed in KS 
cannot be attributed to the sweeping effect. 
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